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I. INTRODUCTION 

Gowdy space-times [l| are solutions of vacuum Einstein equation admitting two com- 
muting, space-like isometries and having closed spatial hypersurfaces. Such models have 
essentially one of the three topologies for the spatial slices: T 3 , S 3 , S 1 x S 2 . The case of 
T 3 is the simplest of these. A further simplification is possible. One can restrict to the 
so-called polarized models in which the two Killing vectors are orthogonal. This simplest 
case of polarized, T 3 vacuum Gowdy model is the focus of this series of works. In this case, 
the complete set of exact solutions is known [2] which generically have initial singularity. 
But there is also an infinite sub-family of solutions for which all curvature invariants are 
finite. The approach to classical singularity is well studied and is known to follow a special 
case of the BKL scenario known as asymptotically velocity term dominated near singularity 
(AVTDS) [3]. At late times, the model is known to be asymptotically homogeneous Q]. 

These models have been analysed in the canonical framework in both metric variables 
as well as in terms of the complex Ashtekar variables. The first attempts of quantization, 
were carried out in ADM variables in jlj] . Another approach which has been more successful 
was based on an interesting property of the model. After a suitable (partial) gauge fixing, 
these models can be described by (modulo a remaining global constraint) a "point particle" 
degree of freedom and by a scalar field </> which is subject to the same equations of mo- 
tion as a massless, rotationally symmetric, free scalar field propagating in a fictitious two 
dimensional expanding torus. This equivalence was used in the quantization carried out in 
6J. Subsequent analysis has been carried out in a large number of works some of which are 
listed in [7]. However in these quantizations, the evolution turned out to be non- unitary 
and in 8| a new parametrization was introduced which implemented unitary evolution in 
quantum theory. 

Canonical description of unpolarised Gowdy T 3 model in terms of the complex Ashtekar 
variables has been given in A complete set of Dirac observables is also known [h]]. The 
canonical quantization of this model was carried out in llj and the physical Hilbert space 
was obtained. Although quantization has been carried out, the difficult issue of singularities 
has not been fully addressed (see however [12J for preliminary attempts). 

In this series of works, we aim to carry out loop quantization of the polarized, T 3 Gowdy 
model, obtained via symmetry reduction in terms of real connection variables. In this paper, 
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we report the first step of recasting the Gowdy model in the (real) connection formulation, 
including the restriction to the polarized model. In subsequent papers the quantization will 
be carried out and issue of singularities will be addressed. 

In section [TT] we discuss, in brief, Gowdy models and then restrict attention to the po- 
larised T 3 case in metric variables. The form of the metric and the space-time solutions 
are discussed. In section IIHI the unpolarised model is described in terms of real Ashtekar 
variables and a consistent reduction is carried out to obtain the polarised model. Here con- 
sistency refers to preservation of the diagonal form of the metric under the Hamiltonian 
evolution. In section IIVI we explain a gauge fixing leading to recovery of the standard form 
of the solutions. Section |V] summarizes the results and includes preliminary comments. 



II. POLARISED GOWDY T 3 MODEL IN METRIC VARIABLES 



Gowdy space-times jl] are globally hyperbolic solutions of the vacuum Einstein's equa- 
tions which are isometric under the action of the Abelian group T 2 which acts on the spatial 
slices assumed to be closed. This means that there are two independent commuting spatial 
Killing vectors. This condition along with the Einstein's equations restricts the allowable 
choices of spatial topology to be only T 3 , S 3 and S 2 x S 1 (or certain manifolds with one of 
these as cover). If, in addition we make an additional assumption that the Killing vector 
fields which generate the T 2 isometry can be chosen to be mutually orthogonal everywhere, 
we get the so-called polarized Gowdy model. The metric, is then diagonal and can be written 
as [2]: 



ds 2 



,2a/ 



-dT z + d6 2 ) + T{e 2W dx 2 + e~ 2W dy 



-2W, 



(1) 



where d/dx and d/dy are the two Killing vectors and a and W are functions of T and 
periodic functions of 6. 

The solutions for W can be obtained from the second order differential equation: 

d 2 W 



d 2 W 1 dW 
+ 







dT 2 T dT d9 2 
Given a solution, W(T, 6), the function a(T,6) has to satisfy: 



(2) 
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The W equation (j2J) encodes the dynamics while the ([3]) encodes the constraints. Inciden- 
tally, this makes the initial value problem and the problem of preservation of constraints in 
numerical relativity trivial for this model. The initial values of the dynamical variable W 
can be freely specified and given a W the constraint a can be trivially determined [3] . The 
requirement of a being a periodic function of 9 imposes a condition of the solutions of the 
W equaW n am ely, J VOrWO.W = 0. The genera, solution gl ven below, does satisfy this 
condition [2|]. 

The general solution to (|2J) is given by 

oo 

W = a + (3lnT + J2 [ a n Jo {nT) sin {n9 + 7n ) + b n N {nT)sm{n8 + S n ) (4) 

71=1 

where a, ft, a n , b n , ^ n and 5 n are real constants and J and N are regular and irregular 
Bessel functions of the zeroth order. The special case of homogeneous model is given by 
ft — | and a n =0 = b n and corresponds to the flat Kasner solution described as: 

ds 2 = -dT 2 + d6 2 + T 2 dx 2 + dy 2 (5) 

It can be shown that the curvature invariant C = R a b c dR abcd blows up almost everywhere 
as T — > 0+. The solutions are therefore generically singular. However for the special choice, 

b n = , ft = \ (6) 

the curvature invariant remains bounded and all components of R a bcd have finite limit as 
T — > 0+. It can also be shown that these nonsingular solutions are analytically extendible 
0] but are causally ill-behaved (have closed time-like curves) in the extended portion. Thus 
there exists a infinite number of nonsingular solutions which however form a set of measure 
zero in the space of solutions. Curvature unboundedness is the generic behaviour. 

III. POLARISED GOWDY T 3 MODEL IN REAL ASHTEKAR VARIABLES 
A. Unpolarised Case 

We begin with the connection formulation in the notation of [ijj], with P l a = E^/^K'j) 



substituted. The classical symmetry reduction to t 
of the complex connection has been studied by 



le unpolarized Gowdy T model in terms 
, ■ Translated into real variables, the 



symmetry reduction is achieved by setting to zero the following components of the densitized 
triad and connection [15I ]: 



E° = = E% , A\ = = A 3 p ; p = x, y ; I = 1, 2 . (7) 

In these variables, the Gauss, the diffeomorphism and the Hamiltonian constraints jl^j ]. 
are given by (k := 87rG Newton , G := G 3 , C := C e ): 

G = i- [deEl + e/A^] ; e/ := e 3 / (8) 

C = — [(^4)^ + ej K A J p E^ - « 7 i4|G? 3 ] 5 (9) 
fey 

^ = ^^=^ [ 2 ^^3<^ + A}E>A$E° K - AfE P jAtE° K - 2ef (d e A J p )E p K E° 

- (1 + 7 2 ) (2KlEtK J p E P j + K J p E P jK«E° K - Kf E p K J a E° K )] . (10) 

In the above, K % a are the components of the extrinsic curvature which are related to the 
gravitational connection A l a and the torsion- free spin connection, r*, as: K l a = / ~f~ 1 (A t a — T l a ) 
and 7 is the Barbero-Immirzi parameter. The spin-connection is defined in eqn. ( 1A.2I) . 

Since none of the quantities depend on x or y we can integrate over the T 2 and write the 
symplectic structure and the total Hamiltonian as: 

n = J d6 (dA 3 e A dE e 3 + dAj A dE p ) (11) 
fltot = 4vr 2 J d6 {X 3 G + N e C + NH) (12) 

Under the 9 coordinate transformation E% transforms as a scalar, E p, s transform as scalar 
densities of weight 1, A z d transforms as a scalar density of weight 1 and A ! p s transform as 
scalars 1 . 

For each p, the A p and E p , rotate among themselves under the U(l) gauge transformations 
generated by the Gauss constraint. It is however possible to choose variables which are gauge 
invariant and will turn out to be more suitable for loop quantization (see section |V]). These 
are introduced through the following definitions: 

El = E x cos(3 ; E x 2 = E x sinf3 (13) 



1 In one dimension, under orientation preserving coordinate transformations, a tensor density of contravari- 
ant rank p, covariant rank q and weight 10, can be thought of as a scalar density of weight = w + q — p. 
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E\ = -E y smp ; E\ = E y cosP (14) 

A\ = A x cos(a + 13) ; A\ = AeSin(a + 0) (15) 
Al = -A y sin(a + (3) ; A 2 y = A y cos(a + (3) (16) 

The angles for the connection components are introduced in a particular fashion for later 
convenience. 

The radial coordinates, E x , E v , A x , A y , are gauge invariant and always strictly positive 
(vanishing radial coordinates correspond to trivial symmetry orbit which is ignored). 
In terms of these variables, the symplectic structure ffTTl) gets expressed as: 

4vr 2 



K7 



7r 2 r 

— J d9 [dAl A dE e 3 + dX A dE x + dY A dP 5 ' + d(3 A dP^ + d/3 A dP' 3 ] (17) 



where: 



X := A x cos(a) ; K := A y cos(a) (18) 

p/3 ;= -E x A x sin(a) ; P^ := -PM y sin(a) (19) 

The gauge transformations generated by the Gauss constraint shift /3, j3 rendering a and a 
gauge invariant. From now on we will absorb the An 2 and use k' := — 2GN ° wton ■ 
It is convenient to make a further canonical transformation: 

£ = ; v = P + P (20) 

p/3 _ p/3 pj3 , pf3 

P" = ^f- i ^=^- ( 21 ) 

In terms of these variables the Gauss and the diffeomorphism constraints can be written 

as: 

G = — \d $ E% + 2P r 'l (22) 
C = — [(d e X)E* + (d d Y)E y - (d 6 E°)Al + («%r/)P" + (3,0 (23) 

The Hamiltonian constraint is complicated but after putting Kf = (Af — Tf)/j, substituting 
the explicit expressions of T", and further simplification, turns out to be: 

r__J_ 

2k Je 

+ 2P"(«%lnP 3 e + (ta<) dg() \ + 2<| (cos 2 (XE x YE y + (P") 2 - (P 5 ) 2 ) 



H= - 



El {(XE X + YE y )d e r] + (XE X - YE y )d B £ - 2P^d e (^ n ^j 
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+ (XE X + YE y )E 6 z A^ + (sm2£){(XE x + YE V )P^ - (XE X - YE^P^ 



cos / V ( cos 

- hzr ) (24) 

where £ = \E° E x E y (cos£)\. 

Under the action of the diffeomorphism constraint X, Y, E°, rj and £ transforms as scalars 
while E x , E y , A z e , P v and P$ transform as scalar densities of weight I. 

This completes the description of the unpolarised Gowdy T 3 Model in the variables we 
have defined. The number of canonical field variables is 10 while there is a 3-fold infinity 
of first class constraints. There are therefore 2 field degrees of freedom. We now need to 
impose two second class constraints such that the number of field degrees of freedom are 
reduced from two to one (as it should be in the polarized case). 



B. Reduction to Polarized model 

The spatial 3-metric, g a b '■= e a e l^ij, with the co-triad e l a defined through e l a E°- := 
5)VE, e\E\ := 6 b a VE , is given by 

ds 2 = C0< ^J_ d6 2 + -^L^-dx 2 + AriUy 2 - 2^-si< dxdy (25) 
s El cos^E x cos^Ev y cos£ ^ y v > 

For the Killing vectors d/dx and d/dy to be orthogonal to each other, the dxdy term in 
the metric should be zero. This implies that the polarization condition is implemented by 
restricting to £ = sub-manifold of the phase space of the unpolarized model. For getting a 
non-degenerate symplectic structure, one needs to have one more condition. This condition 
should be chosen consistently in the following sense. 

We expect the two conditions to reduce a field degree of freedom. This can be viewed 
in two equivalent ways. The condition £ = makes the metric diagonal and this property 
should be preserved under evolution (i.e. the extrinsic curvature should also be diagonal). 
Alternatively, the unpolarized model is a constrained system and we want to impose two 
conditions such that one physical (field) degree of freedom is reduced. The extra conditions to 
be imposed should therefore be first class with respect to the constraints of the unpolarized 
model i.e. should weakly Poisson-commute with them. 
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Indeed, this can be done systematically by viewing £ = as a new constraint 2 and 
demanding its preservation under the evolution generated by the total Hamiltonian. Since 
£ = weakly Poisson commutes with the Gauss and the diffeomorphism constraints, only 
the Poisson bracket with the Hamiltonian constraint is needed. 

m « , {£(6), J dd'N{e')H{d')} « . (26) 

It follows that, 

£(6)^0 x(0) ■■^2P^ + E d 3 d e (]nE y /E x ) wO (27) 

The Poisson Bracket of \ with the Hamiltonian turns out to be zero on the constraint surface 
i.e. x ~ X ~ 0. Thus, the reduction to the Polarized model is obtained by imposing the 

two polarization constraints 

^0 ; x~0 ; {m,x(0')} = 2^5(9-9') (28) 

Remark: To see that the x ~ condition follows from preservation of g xy = 0, note that, in 
the metric formulation, for the present case, it implies that g xy ~ K xy = K x e y {= e x K y ) = 0. 
Using the definition K l a = , ~f~ 1 (A , a — T^) and the expressions given in the appendix, one 
can check directly that K xy = x = 0. Note that this is not equivalent to requiring 
orthogonality of components of the connection, A l x A % = which would imply a = a (see 



eqns. (|T5l and [TBI) . This condition, mentioned in the literature [15|, [16|], is very different from 
the x ~ condition and is not preserved under evolution. 

It follows from (|22l) that {x, G} = And using (|23|) . one can see that : 

j^yVc^} = N%^0 , jx^iV^} = d e (N d X )^0. (29) 

We can solve the polarization constraints strongly and use Dirac brackets. Symbolically, 

{/, qY = {/, ^} - {/, ei © xr 1 © {x, ^} - {/, x} © {x, e} -1 © {£, ^} 

Here denotes appropriate integrations since we have field degrees of freedom. 

Since the polarization constraints weakly commute with all the other constraints, the 
constraint algebra in terms of Dirac brackets is same as that in terms of the Poisson brackets 



2 The choice £ = also requires E$ > for the spatial metric to have signature (+, +, +). The choice £ = it 
would require E\ < 0. From now on E\ > will be assumed. 



8 



and thus remains unaffected. Furthermore, equations of motions for all the variables other 
than £, P^ also remain unaffected. We can thus set the polarization constraints strongly 
equal to zero in all the expressions and continue to use the original Poisson brackets. 

The expressions of the constraints simplify greatly and in particular the Hamiltonian 
constraint simplifies to, 
T 2 1 



H = 



2k y/E 



\2 



2 



(XE X + YE y )E 9 3 d e r] + 2^XE x YE y + (XE X + YE y )E°A 3 e } 



2 



{{d e ElY-(EldeHEy/E x )) 2 ) 



where, we have also eliminated P v in terms of the Gauss constraint using 2P V = (K,^G—d e E 3 ) 
and E = \E 9 3 \E x E y . 

Noting that 77 is translated under a gauge transformation, we can set it to any constant 
and fix the gauge transformation freedom. Explicitly, imposing 77 m as a constraint, we 
can fix the A 3 from preservation of this gauge fixing condition. Once again we can use Dirac 
brackets with respect to the Gauss constraint and the 77 ~ constraint and impose these 
constraints strongly. With this done, the first two terms and the G dependent piece in the 
last term in the Hamiltonian, drop out and so do the degrees of freedom 77, P v . We are left 
with six canonical degrees of freedom and the two first class constraints, leaving one field 
degree of freedom. Thus our final variables and constraints for the polarized Gowdy model 
are (absorbing away the Immirzi parameter): 

« := ,£:=E°,A := 7^ , K x := 7^ , K v := (31) 

{K x (6), E x (6')} = k5(6 - 6') , (and similarly for (K y , E y ), (A, S) pairs); (32) 

C = - [(d e K x )E x + (d e K y )E y - (d e S)A] (33) 



K 



1 



^{(K X E X K V E») + (K X E X + K y E y )£A 
±={(d e £) 2 -(£de(HE y /E x ))) 2 } + d t 



Ed E 



Ay/E^ J V. \[E 

The constraint algebra among the C[N e ], H[N] is: 



(34) 



{ C[N e ] , C[M e ] } = C [ N e d e M e - M e d e N e ] (35) 
{ C[N e ] ,H[N] } = H[ N e d e N } (36) 
{ H[M] , H[N] } = C [ (Md e N - Nd e M)E 2 E' 1 ] (37) 



Since each term in the Hamiltonian constraint is a scalar density of weight +1 and each term 
in the diffeomorphism constraint is of density weight +2, the first two brackets are easily 
verified. The last one also follows with a bit longer computation. We have thus verified the 
constraint algebra of polarized model showing the consistency of the reduction procedure. 

IV. SPACE-TIME CONSTRUCTION 

The next task is to find the set of gauge inequivalent solutions of the Hamilton's equations 
of motion, satisfying the two sets of constraints and obtain the space-time interpretation. 
The total Hamiltonian being a constraint, the Lagrange multipliers - the lapse function and 
the shift vector - also enter in the Hamilton's equations of motion. These need to be either 
•prescribed or deduced via a gauge-fixing procedure. Once this is done, one can obtain the 
solution curves in the phase space with "initial points" lying on the constrained surface. 
The space-time metric, solving Einstein equation is then given by, 

ds 2 = - N 2 (t, x { )dt 2 + gij (t, x 1 ) (dx { - N*(t, x { )dt) (dx j - N j (t, x^dt) (38) 

For our case, the metric is diagonal, (x 1 , x 2 , x 3 ) <-> (9, x, y), N" 1 (iV e ,0, 0) and the metric 
is independent of the coordinates (x,y). The t = constant, hyper- surfaces are diffeomorphic 
to the 3-torus. The metric components are given by: gee = E x E y £~ 1 = E£~ 2 ,g xx = 
SE y /E x ,g yy = SE X /E y (eq. (125]) ). The Gowdy form of the metric ([T|) is realized if one 
prescribes N e = and N 2 = gee- 

Such a prescription is eminently consistent since any metric on a two dimensional man- 
ifold (coordinatized by (t,9)), can always be (locally) chosen to be conformally flat. This 
however does not fix the coordinates t, 9 completely - one can still make the conformal dif- 
feomorphisms: t —>■ t' = t + 9), 9 — > 9' = 9 + £, e (t,9) with £ satisfying the conformal 
Killing equations: d t ^ — dei 8 = = — d t £ e . 

We will first take the above prescription for the lapse and the shift, obtain the Hamilton's 
equations of motion, use the freedom of conformal diffeomorphisms and reduce the equations 
to those given in section II. Subsequently, we will also exhibit gauge fixing functions to 
arrive at the same result. This will complete the identification of inequivalent solutions of 
the Einstein equation. 
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With the choices N e = 0, N = y/EE' 1 , the space-time metric (1381) is 

f E y E x \ 

ds 2 = E£~^-dt 2 + d9 2 )+£(=-dx 2 + =-dy 2 ) , (39) 



E x Ev 

and the time evolution is governed by the Hamiltonian alone which is given by, 
H[£- x Ve] = \f de -^{(K x E x K y E y ) + (K x E x + K y E y )£A^ 

VE n f£d e £ 



1{(^) 2 - (£d e (HEy/E*))f } + *fde 



A£V ° > v - ^ ' J £ V ^/~E 

In anticipation let us define 2W := ln(E y /E x ) and 2a := hi(E x E y /£). One obtains 



(40) 



J = £-\K y E y + A£) , ^- y =£-\K x E x + A£) , £ = (K X E X + K y E y ) (41) 

2d t W := ft ln^ = (j^gB (42) 
E x E y 

2d t a := d t In— = 2 A (43) 
The Poisson brackets of K X E X , K y E y with the Hamiltonian are given by, 

{K X E X , H[£^VE}} = l -d e (^ln|^ + l -d 2 e £ (44) 

{K y E y , H[£-*VE\} = --d 6 (Wn— J + -d 2 £ (45) 

{i^ - K y E y , H\£- 1 Ve}} = 8 e Ud e \n^\ (46) 
{K X E X + K y E y 1 H[£- l VE}} = d 2 e £ (47) 

From these, we get second order equations for £ , W as, 

d 2 £ = d 2 £ (48) 
d 2 W = ^d e (£d e W) - [k^j 9 t W (49) 

The equation for £ is a simple wave equation and given a solution of this, the equation 
for W can be solved determining W or the ratio E y /E x . From the first order equations, 
one determines the K X E X ± K y E y as well. The Hamiltonian constraint then determines A 
in terms of known quantities, £, W and the ^-derivatives of a. Using the equation d t a = A, 
one obtains, 

9 <- A -If + h m + {mf) ~ % d ' a ~ I w + % (50) 
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One can also obtain, by direct computation and using the diffeomorphism constraint, 



d t £ 



2d t Wd e W - -^-&a + 



d e £ Q , d%£ d t Sd e S 



(51) 



£ 1 £ 2£ 2 

From these two equations, one can obtain dga, d t a in terms of £ and W which can be 
integrated. Thus the metric can be completely determined starting from a solution for £. 

However, all these solutions are not gauge inequivalent corresponding to the fact that the 
coordinates can still be subjected to conformal diffeomorphisms. Under these coordinate 
transformations, £ which is the determinant of the metric on the symmetry torus, is a 
scalar. Under conformal diffeomorphisms, the wave operator gets scaled by a prefactor. 
Hence, under the transformations generated by conformal Killing vectors, solutions of the 
wave equation transform among themselves. In fact the conformal Killing vectors also 
satisfy the wave equation and on the cylinder (t, 8), both £ and £ satisfy the same boundary 
conditions. Thus, their general solutions are linear combinations of exp{in(t ± 6)},n ^ 
and a solution of the form A + Bt. The Killing vectors however satisfy first order coupled 
equations. This removes the ^-independent, linear in t piece from the general solution. 
Consequently, one can use conformal diffeomorphisms to remove the ^-dependence from the 
solutions for £ as well as the constant piece. In other words, all solutions for £ , except 
£ = #t are related to each other by conformal diffeomorphisms. The gauge inequivalent 
solutions are thus obtained from the choice £ = t. Equivalently, one has finally fixed the 
(t, 6) coordinates completely. The time coordinate so fixed will be denoted by T. 

With this choice, £ = T, the constraints also simplify to, 

= £(d e W) 2 - £- x {K x E x K y E y + (K X E X + K y E y )A£} (52) 
= E x d e K x + E y dgK y , (53) 

one gets K X E X + K y E v = 1 and the equations f|49l [50| IBTj) go over to the equations ([21 E]) . 
From these one recovers the usual solutions listed in section II. 

Remark: Up to the derivation of the equations for £ and W, the constraints are not 
used. The T 3 topology has also not been used! Thus these expressions are also valid for 
polarized versions of Gowdy models with other topologies. In Gowdy's original analysis, the 
three allowed topologies are distinguished by different choices of solutions of the equation 
for £ (R, the determinant of the two metric on the T 2 orbits, in Gowdy's notation). The 
different topologies get distinguished by the boundary conditions on £ and on the conformal 
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Killing vectors. For non-T 3 topologies, 9 G [0, n] and have to vanish at the end-points. 
With these taken into account, the gauge inequivalent solutions are obtained by choosing 
£ = sin(t) sin(0) fl. 

We reproduced the known results by obtaining the solutions of the Hamilton's equations 
with chosen lapse and shift, motivated by comparison with the space-time form of Gowdy 
model, and invoking the 'residual' freedom in the space-time coordinates to obtain the gauge 
inequivalent solutions. Thus we used the canonical structure as well as the anticipated form 
of space-time geometry to arrive at the distinct solutions. We would like to see if the same 
result can also be derived by using only the phase space view. 

Within a phase space view, the lapse and the shift are to be determined by doing an 
explicit gauge fixing. To do this, we will keep the lapse and the shift as unspecified and look 
at the evolution generated by the total Hamiltonian, 

H tot [N 9 , N] = d9N e {E x d e K x + E»d e K y - Ad S} 

(K x E x K y E y ) + [K X E X + K y E v )£A 



al9N 



^{(deS) 2 - (Sde(HE y /E x ))f } + d e 



(54) 



4a/E I J \\[E 

Denoting by over-dots, the Poisson brackets with the total Hamiltonian, it is straight 



forward to see, 

™ , 9,{N»E') 



E x 


N 


E* 




& 


N 


E~y 




£ 


N 


£ 


Ve 


(K X E X ) 


= -de 


(K y Ev) 


= l -de 
2 



Ev 
£ 



(55) 
(56) 



(K X E X + K y Ey) + — 5i (57) 



^| (edging + deS^j } + d e (N K x E x ) (58) 
{-^= (-Sde^ + de£) } + d e (N e K y Ey) (59) 



The following combinations are convenient for looking at gauge fixing. 



(K X E* + K y Ev) = d e {^<^} + d e { N ° ( K ^ X + K y E y )} (60) 
N£ 

£ = -= (K X E X + K y Ey) + N 9 d e £ (61) 
v E 
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f N£ E y 1 

(K X E* - K y Ev) = d e l—£d e \n—\ + d e {N 9 (K X E X - K y E»)} (62) 
/ w.y\ N E y 

(in-) = — (K X E X - K yE y) + N%ln- (63) 

The first two equations above show that we can consistently impose K X E X + K y E y = Ci, 
a constant, and d$£ = as two gauge fixing conditions. Preservation of the first leads to 
N 9 = f(t) while that of the second leads to N£/\^E = g(t). Since d$£ = already requires £ 
to be a function of t alone, we can strengthen the gauge fixing condition by specifying £ = t. 
This determines N = C\\/~E£~ X . Evidently, we must have a non-zero lapse and therefore 
Ci 7^ must be chosen. The sign of C\ will determine if £ increases or decreases with t 
and by convention we can take the sign to be positive and without any loss of generality, we 
choose Ci — +1 and denote the t by T as before. 

The shift is however determined to be a function of T alone. With such a shift, CfiV 6 '] = 
f(T) J C generates T-dependent translations of the ^-coordinate. All tensor densities on the 
spatial slice, transform as scalars under these translations, and there is no way to fix the 
left over constraint f C, by any gauge fixing condition. However, we can always redefine 
the ^-coordinate such that al6 — f(t)dt =: d6' . This means that solutions inequivalent with 
respect to translations, can be determined by effectively choosing shift = 0. Incidentally, 
for other admissible topologies, the shift has to vanish at 9 = 0, 7r and hence f(t) — is the 
only admissible solution. We have thus achieved our goal of determining the same lapse and 
shift, by explicit gauge fixing. The inequivalent solutions are then obtained as in section II. 

One can make the physical degrees of freedom explicit by noting that 2W = \n(E y / E x ) 
and ir w := K X E X — K y E v are canonically conjugate. Similarly, 2a := — \n(E x E v ) and ir a := 
K X E X + K y E v are also conjugate variables. The gauge fixing conditions are: £ = T, K a = 1 
while the gauge-fixed form of constraints become: 

C = - [ir w d e W + d e d] (64) 

K 



{t-*Ve)h = \ 



AT 



A + T(d e Wy 



(65) 



The Hamiltonian constraint determines A completely in terms of W, t\w while the diffeo- 
morphism constraints determines the d except for the homogeneous (^-independent) part. 
The periodicity of a also requires the J ir w d e W = which is a constraint on the W, n w . 
The physical degrees of freedom are thus described by W, ir w together with one constraint 
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and the homogeneous pieces of a, 7Tq. Our gauge fixin g ha s fixed the homogeneous part of 



Hi 



TT a to be 1. These are of course the well known results 

Observe that in the homogeneous limit (all variables independent of 9), one gets the 
Bianchi I model. The Hamiltonian constraint, for each 9 looks like a Bianchi model with a 
potential and is highly suggestive of the BKL scenario and has been explored numerically 
as well [17 ]. 



This completes the canonical formulation of the polarized Gowdy model on T 3 in terms 
of the real connection variables. 



V. DISCUSSION 



In this paper, two main reformulations of the polarized Gowdy model in real con- 
nection variables have been done. First is the choice of the gauge invariant vari- 
ables: A x , A y , E x , E v , a, at and the subsequent canonical transformation to the variables 
X, Y, P^,P V . This has already been done in the case of spherical symmetry and also men- 
tioned for cylindrical waves in [lj|. The main advantages of these variables are that the 
volume becomes a functional of the momenta variables alone and the components of the 
connection along the homogeneous directions are separated neatly and gauge invariantly, 
into extrinsic curvature components (X, Y) and the spin-connection components (T x ,r y ). 
In the quantum theory, both the features allow a simpler choice of edge and point holonomies, 
simpler form for the volume operator and also a more tractable form of the Hamiltonian 
constraint fis| . 

The second aspect, obtains the polarized model from the unpolarized by a simple sys- 
tematic reduction (Dirac procedure) ensuring a consistent reduction at the level of physical 
degrees of freedom. Getting this reduction consistently is important since the form of the re- 
duced constraints, depend on the reducing conditions. In contrast to the second polarization 
condition mentioned in the literature, namely orthogonality of the connection components 
in analogy with that of the triad components, our x ~ condition, (|27j) is consistent with 
dynamics. The consistency is seen in three ways: from a systematic derivation, verifying 
the constraint algebra of the reduced constraints and finally reproducing the known space- 
times, obtained by directly solving the Einstein equations for polarized ansatz. We are thus 
confident to use these constraint expressions in the passage to quantization. 
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We would also like to note that in the reduction to polarized model, we had two options: 
£ = {El > 0) or £ = 7r {El < 0). In the metric variables and classically, either one 
of these suffices. (In the triad variables, these two correspond to opposite orientations.) 
The subsequent gauge fixing was also naturally restricted to one of these choices (we chose 
the former). At this stage, one could imagine doing a "loop quantization" of the gauge- 
fixed model which now has a true Hamiltonian and explore the fate of the singularity. In a 
quantum theory however, one could have an extension across the degenerate triad and this 
will be missed in a quantization of the gauge-fixed model. 
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APPENDIX 

In this appendix, we collect some of the useful expressions to help reproduce the compu- 
tations in the main text. 

The symmetry reduction leaves only the following non-zero components of the gravita- 
tional connection and the densitized triad: 

A\ — > AlA 1 ^ 1=1,2 p = x,y , 

E? — E e 3 ,E p j, 1=1,2 p = x,y. (A.l) 

The triad has the components: e? = E^E^ 1 ^ 2 with E := detE? = E 6 Z A, A := EfE%-E%E%. 
The co-triad (inverse triad), e l a , has the components: e\ = y/A/El, el = y^El/A E\, e* = 
-^Ej/A El, e\ = -^Ef/A E\, e\ = ^Ef/A Ef. 
The spin connection is defined by, 

:= - e^ej [d [a e k h] + \e c k e l a d [c e l b }j . (A.2) 

Of these, T 3 p = = are identically zero. The remaining components are given by, 

11 = T^r (EfdeEf - E\d e E\ + E*d e E y 2 - E y 2 d e E x 2 ) 



2A 
1 



r 1 = . 

2V A 





E e \ E e - - / E x \ E e - - / E y \ 

-t E " ) + ^ E "- Evd » (ti Evd « (ti) 
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r 1 

V 




' Ell 1 - -2-Ev ■ Evd e 



E 

- -JlE x ■ Evd e 
A 

E 9 

+ -^E x ■ Evd e 



El 



El 



E 
E\ 



E 

+ -^E x ■ Evd e 



El - - 
+ 3 ^ . Ex q 

A 

E e 

- -?-E x ■ E x d e 
A 



E\ 



(A.3) 



E\ 



E 
E\ 



E%E% etc. 



where, E x ■ Ev := EfE\ 

In terms of the radial and angular variables E x , E y , E®(= £),£,rj given in equations (fT3l 
HEOD, one has A = E^cosf , E x ■ E x = (E x ) 2 , & ■ & = (Ey) 2 , E x ■ & = EFE»wn£. 

In the computation of the Hamiltonian constraint, one needs the combinations: EfY\ 
and E^A l b . The non-zero ones are given by, 



EX 

E?ri 



EfAl 



E X A\ 



EfAl 



-E% 



E y 

tan£ d e ( In— ) - d e r] 



E x 



-d e (£ 3 e tan£) 



~\ de 
K X E X 



E x El 
Ev cos£ 



E y r 

i y 



E y r 



- -do (£ 3 e ta<) 



5* 



KyE y 



E y El 

E x cos£ 
E?At = 



(A.4) 



(jpj {-(J* + P")cos£ + K x E x sm£} 

(^) {(p"-p«)co S e+^si<} 



(A.5) 



The Hamiltonian constraint in eqn (flOl) is simplified by eliminating the extrinsic curvature 
in terms of the gravitational connection and the spin-connection, K\ := 7 _1 (t4^ — r*), and 
using the above equations. One begins with the expression: 
1 1 



H 



2k J~E 



1 1 

2k ^E 



[ e ijk E«E){d a A$ - d b A k a ) + e ijk e%E? E)A™A n b (A.6) 
-(1 + 7 " 2 ) {E?(Ai - r a )E){Ai - T{) - Et{A{ - ^^(4 - P 6 )}] 
[ {e ljk EtE]{d a A k b -d b A k a )} (A.7) 



-7" 2 {{E«A\){E]Ai) - {E?Al){E]Ai)} 
-(1 + T 2 ) {{E«r a ){E)T{) - (Etr b )(E b fi)} 
-2(1 + 7 - 2 ) {-{EtA\){E]T{) + {EtT\){E]Ai)}] 

The terms quadratic in A's combine to get 7~ 2 , while terms linear in A and T get a prefactor 
of 2(1 + 7~ 2 ). In terms of the angular and radial variables given in ()13p . the terms in the 
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braces in eqn (1A.7j) become: 

First = 2e JK E°E p jd e A K 



Second 



Third 



El [ (K X E X + K y E y )d dV + (K X E X - K y E y )d £ - AdgF" 

+2P r 'd e \n(E x E y ) - 2P^d e ln{E y /E x )] (A.8) 
{ElAl + K X E X + KyEvf - {E 9 3 Al) 2 - {E^A I a ){EjA J p ) 
2 cos 2 (0 {K x E x K y E y + (P") 2 - (P«) 2 ) + 2 P 3 e ^ (K X E X + 
+ 2 sin(Ocos(0 {P^(K X E X + - P r \K x E x - K y E y )) (A.9) 

- (E^l) (E°jT J p ) 



(deEiy 



p 



cos(0 



- 



p 



cos(0 



d ]n{E u /E K 



(A.IO) 



Fourth = -(P 3 ^)(P 3 ^ + K X E X + K y E y ) + (P 3 ^)(P 3 ^) + (P^)(P^) 

= [-2P 6 d e \n(E y /E x ) + 2P r < d \nE 9 3 + 2P r? tan(0<%£ 
+ (K X E X + K y E y )d e n + (P X P* - P.P^H^] 



(A.ll) 



Using these equations and with further rearrangement of terms, leads finally to the Hamil- 
tonian constraint of eqn. (12^1) . 

The polarization conditions are f = = 2P € + P^ ln(P^/P x ). With £ = (3 - (3 = 



= — sm/3 , r 2 = cos/? r a 



ipf 

2E X ' 



E x J 



1 (9 


= -(1/2)^, 






r 1 


= cos/? r y , 


r 2 

2/ 






:= ifftln ( 
2P^ \ 




pa 
p^ 



(A.13) 



Using the Gauss constraint and the polarization constraint, x — 0, it follows that, 



pP pP 

A x sin a := — — — = T x , A v sin a := — — = — r„ . 
E x Ev 



(A.14) 



Notice that T x , T y are gauge invariant. 

In checking preservation of various constraints as well as verifying the constraint algebra, 
the following Poisson bracket is useful, 



K X E X {9) 



d e iE x {6') 
E x (6') 



2G 



k :- 



Newton 



7T 



(A.15) 
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